Optics Letters

The spatial resolution of optical imaging systems is fundamentally restricted by the diffraction limit. However, in widefield live-cell microscopy, the achievable resolution is further constrained by the specimen motion, which indicates the existence of a fundamental spatio-temporal resolution trade-off between signal accumulation during the full frame integration and the resulting motion blur. To improve the fidelity with which moving objects can be imaged, a quantitative understanding of this spatio-temporal trade-off is necessary. Here, we present a systematic analysis of motion-induced resolution dynamics measured with spectral signal-to-noise ratio (SSNR). We developed a simulation framework which models the image formation of objects undergoing arbitrary motion, to evaluate the degradation of the spatial resolution under translational and rotational dynamics. Our results demonstrate that for translating objects, the spatial resolution is anisotropically reduced as a function of the orientation of the object relative to the motion vector, leading to the spectral signal-to-noise ratio degrading by up to 50% and the resolution by up to 40% for a 90{degrees} change in the motion direction. Furthermore, we show that for rotational motion, conventional radially averaged metrics such as the Fourier Ring Correlation are not able to quantify the effects of angular blur. On the other hand, the SSNR is able to accurately quantify this degradation. These findings underscore the necessity of an object-oriented imaging approach, in which acquisition parameters such as exposure time are tuned to specific biological spatio-temporal characteristics to optimize the trade-off between motion blur and spatial fidelity.

Microbubble contrast-enhanced ultrasound (CEUS) relies on discriminating nonlinear bubble signals from linear tissue backscattering. While Singular Value Decomposition (SVD) filtering improves this discrimination, existing techniques often fail to retain the slowly-moving microbubble signals from static clutter. Here, we present a novel multi-stage singular value decomposition (MS-SVD) framework for ultrafast CEUS imaging. Our method employs plane-wave transmissions at multiple angles and acoustic pressure levels (implemented via duty-cycle modulation) and alternating transmit polarity. The beamformed data are then processed by three sequential SVD filters: (1) spatial-angular SVD to extract coherent signals across all transmit angles, (2) spatial-pressure SVD to separate linear fundamental and nonlinear harmonic components, and (3) spatiotemporal SVD to isolate moving microbubble echoes from tissue clutter. In in vitro flow phantoms and in vivo rat brain through a cranial window, MS-SVD dramatically improves microbubble detection compared to conventional SVD filtering, MS-SVD yields much stronger vascular contrast and suppresses tissue clutter to a greater extent. The resulting power-Doppler and super-resolution maps are notably cleaner and more complete: MS-SVD detects substantially more microbubble events in ULM, revealing finer vessel details and more accurate flow speeds. By capturing the full acoustic signature of microbubbles (both fundamental and harmonic), MS-SVD achieves higher contrast-to-noise and sensitivity in CEUS. These gains make it a powerful front-end for super-resolution ultrasound localization microscopy and other high-sensitivity microvascular imaging applications.

Three-dimensional photoacoustic tomography (3D-PAT) enables noninvasive structural and functional imaging with optical absorption contrast and ultrasonic detection depth. However, its spatial resolution is limited by acoustic diffraction, and incomplete detection geometry can substantially degrade image fidelity and quantitative accuracy. Here, we present a ULM-guided model-based reconstruction framework, termed 3D-PAULMprior that incorporates sub-diffraction vascular priors from concurrent ultrasound localization microscopy (ULM) into 3D photoacoustic reconstruction. The method uses weighted regional Laplacian regularization to integrate high-resolution vascular information into the inverse problem, thereby enhancing vascular sharpness, suppressing limited-view artifacts, and improving blood oxygen saturation estimation. We validated 3D-PAULMprior using numerical simulations, tissue-mimicking phantoms, and in vivo mouse brain imaging. Compared with conventional reconstruction, 3D- PAULMprior improved spatial resolution by over 50%, increased contrast-to-noise ratio by 261.2%, and enhanced structural similarity index by 24.6%. In vivo, 3D-PAULMprior recovered vascular structures that were poorly resolved or missing in conventional reconstructions and produced more spatially confined sO2 maps. These results establish 3D-PAULMprior as a robust multimodal reconstruction strategy for high-resolution structural and functional photoacoustic imaging.

For decades, the photon counting histogram (PCH) was used as the sole method to quantify fluorophore numbers in a diffraction-limited focal volume. This technique combines spatial excitation profiles, and the distribution of photon counts to register the photon emission statistics of individual fluorophores. However, this approach has not yet been transferred to widefield fluorescent imaging due to the lack of fast and single photon sensitive camera sensors which can capture the photon emission statistics of a single fluorophore. Here, we explore avenues towards quantitative analysis of the active fluorophore number by leveraging recent advancements in single photon avalanche diode (SPAD) array technology. Binary exposures of a SPAD array can be synchronized with picosecond laser pulses to measure the PCH in a widefield setting. Then, by modeling the statistical relationship between the active fluorophore number and the PCH in a region of interest following a laser pulse, we can perform Bayesian inference of this number. The model is demonstrated experimentally by counting quantum dots and various numbers of fluorescent dye molecules bound to DNA origamis. We find that this method has several important applications in widefield microscopy, including enhanced localization microscopy and constrained fitting of multiple unresolvable fluorescent emitters.

We introduce a workflow to identify oligomeric structures that are recorded with single-molecule localization microscopy (SMLM) under cryogenic conditions. Typically, these oligomers are assumed to consist of protomers arranged as equilateral two-dimensional polygons and every protomer is labeled with a dye molecule for visualization. Unlike previous work, we consider scenarios in which the sample plane has an unknown orientation relative to the focal plane. Our contribution is a high-precision plane-fitting algorithm to determine the sample plane, combined with geometrical transformations and two circle-fitting algorithms to identify the oligomeric structures. Our simulations on synthetic data demonstrate that the proposed workflow achieves high accuracy in estimating both the unknown tilted plane and the oligomer size.

High bit-depth medical images preserve subtle intensity variations that are important for quantitative analysis and clinical interpretation, but their large dynamic range poses challenges for efficient compression. We propose a bit-plane-aware dual-stream compression framework for 16-bit medical images by separately modeling the most significant bit (MSB) and least significant bit (LSB) components. The MSB structural stream is encoded using JPEG coding with a Duplicate Segment Skipping (DSS) strategy to exploit spatial and segment-level redundancy, while the LSB detail stream is compressed using learned image compression to represent residual variations and fine-grained details. Experiments on four MRI and CT datasets show that the proposed method consistently outperforms representative traditional and learning-based codecs, achieving the lowest bit rate across all datasets. Meanwhile, it preserves high reconstruction fidelity. As a downstream application, we further demonstrate that the compressed bitstreams can be effectively integrated with DNA encoding and converted into sequences with favorable biochemical properties.

Stimulated emission depletion (STED) microscopy is a super-resolution fluorescence imaging technique that achieves high spatial and temporal resolution by exploiting stimulated emission to induce fluorescence depletion (FD) and is expected to have substantial utility for imaging applications using fluorescent proteins. However, the compatibility of fluorescent proteins with STED microscopy systems has been understood primarily through empirical observations, and there is no established methodology for the rational selection of fluorescent proteins for STED microscopy. In this study, we systematically evaluated the compatibility of commonly used fluorescent proteins with STED microscopy systems by measuring FD properties using transient absorption spectroscopy and fluorescence dip spectroscopy, both of which are classified as two-color spectroscopy (TCS). Fluorescent proteins identified as compatible with the STED microscopy system based on the TCS measurements were employed for three-dimensional STED imaging of cellular samples expressing each protein. In all samples, three-dimensional spatial resolution was improved relative to confocal laser microscopy, with particularly marked improvements in z-axis resolution. These findings demonstrate that measurements of FD properties via TCS provide a robust approach for evaluating the compatibility of fluorescent proteins with the STED microscopy system and for selecting suitable fluorescent proteins for STED imaging.

Whole-body SPECT bone scintigraphy reflects skeletal metabolic activity throughout the body and plays an indispensable role in the screening, treatment evaluation, and prognostic assessment of bone metastases in tumors. However, the automatic detection and segmentation of hypermetabolic bone lesions remain challenging due to low contrast, limited spatial resolution, and complex lesion distributions. In this study, we proposed Bone-Segnet, a dual-view guided automatic segmentation network for hypermetabolic bone lesions that integrated multi-scale feature modeling, global context modeling, and view-conditioned modulation. Pixel-level annotated anterior and posterior whole-body bone scintigraphy images were used for model training and prediction. The proposed network enhanced the recognition of low-contrast and small-scale lesions through small-lesion enhancement and multi-scale contextual modeling. A Transformer module was further introduced to strengthen global feature representation, while cross-view collaborative modeling was achieved by incorporating the complementary characteristics of anterior and posterior imaging. Experimental results demonstrated that the proposed method outperformed existing approaches across multiple evaluation metrics, with the Dice score improving from 0.7440 to 0.8750, indicating a substantial improvement in segmentation performance. Further quantitative analysis based on the segmentation results revealed significant differences among disease types in lesion count, pixel burden, and spatial distribution patterns, reflecting the heterogeneity of disease-related skeletal metabolic activity. Overall, the proposed method improved automatic lesion segmentation performance and enabled quantitative analysis of lesion burden and spatial distribution patterns, providing objective data support for the assessment of related diseases. Index Terms--Whole-body SPECT, bone lesion segmentation, dual-view modeling, quantitative analysis.

Fast extraction of physically relevant information from images using deep neural networks has led to significant advances in fluorescence microscopy and its application to the study of biological systems. For example, the application of deep networks for kernel density (KD) estimation in single-molecule localization microscopy (SMLM) has accelerated super-resolution imaging of densely labeled structures in the cell. However, localization of fluorescent molecules in dense images is a difficult inverse problem with potentially multiple solutions. To model a probability distribution of solutions to this problem, we propose a generative modeling framework for KD estimation in SMLM based on a conditional variational diffusion model (CVDM). In this framework, CVDM is trained to perform localization tasks on low-resolution measurements by modeling a distribution of high-resolution KD estimates. This approach allows us to probe the structure of the distribution on KD estimates and express uncertainty, which is not currently offered by existing deep models for localization microscopy. We demonstrate that this model permits high-fidelity super-resolution, enables the uncertainty estimation of regressed KD estimates, and has important implications for image restoration in single-molecule and super resolution microscopy.

Background: Quantitative lipid assessment is central to identifying rupture-prone coronary plaques and represents a therapeutic target for lipid-lowering therapy. Near-infrared spectroscopy (NIRS)-derived lipid core burden index (LCBI) is well validated and widely used for detecting lipid-rich lesions. Optical frequency domain imaging (OFDI) is increasingly adopted for guiding percutaneous coronary intervention (PCI) due to its high-resolution structural imaging capabilities. Depolarization-sensitive OFDI (depOFDI) provides intrinsic lipid contrast and may enable combined structural and compositional plaque characterization within a single OFDI-based platform. Objective: To define an OFDI-derived lipid metric and evaluate its agreement with NIRS-derived LCBI. Methods: Thirty-three patients underwent both polarization-sensitive OFDI and NIRS-intravascular ultrasound imaging during PCI. After exclusion of 4 datasets, 29 co-registered pullbacks were analyzed. A signal-to-noise-corrected depolarization metric was used to identify lipid-rich regions and generate depOFDI chemograms. maxLCBI4mm value and location, as well as total LCBI, were computed and compared with NIRS. Results: depOFDI demonstrated strong agreement with NIRS, showing high correlation for maxLCBI4mm (r^2 = 0.862) and total LCBI (r^2 = 0.867), along with strong spatial concordance for the location of the maxLCBI4mm (r^2 = 0.900). Bland-Altman analysis of LCBI4mm showed minimal bias (10.7) with 95% limits of agreement of [81.4 to 102.8]. Conclusions: depOFDI enables accurate quantification of lipid burden alongside the high-resolution structural information inherently provided by OFDI. Because depolarization metrics can be derived from polarization-diverse detection available in many commercial OFDI systems, this approach provides a practical pathway toward comprehensive plaque characterization within existing PCI workflows, without the need for additional imaging modalities.

In previous studies by the author on binocular vision with the asymmetric eye (AE), which models a healthy human eye with misaligned optical components, the results were primarily presented in the Rodrigues vector (RV) framework and supported by simulations and 3D visualizations in GeoGebras dynamic geometry environment. In this paper, the novel geometric kinematics of the human eye, that is, the eye with misaligned optics, and simplified assumptions about the eye rotations (the eyes translational movements are disregarded), are developed within the framework of rigid-body rotations. The originality of the analysis lies in a precise geometric decomposition of a full rotation of the eyes posture into a torsion-free rotation (the geodesic part) and a torsional rotation (the non-geodesic extension of the geodesic part). This decomposition is extended to the corresponding decomposition of the angular velocity. A novel derivation of the eyes angular velocity from the RV formulation of the eye kinematics is proposed.

Visual prostheses face a critical miniaturisation challenge: converting photoreceptor signals to biologically appropriate retinal ganglion cell (RGC) stimulation patterns within the spatial constraints of intraocular implants. Existing systems rely on external microcontrollers for signal processing, limiting scalability for high-density pixel arrays. This paper presents an integrated per-pixel circuit architecture that directly converts photocurrent into frequency-modulated current pulses that match RGC activation thresholds. The design targets are established through NEURON computational modelling of red-green colour-opponent midget RGCs, identifying stimulation thresholds of +0.1nA to +3.5nA for depolarisation and -0.1nA for repolarisation. The proposed circuit combines a transimpedance amplifier, a voltage-controlled oscillator with a Schmitt trigger, and a current-controlled output stage to generate biphasic pulses within these thresholds. A complementary output provides lateral inhibition, reducing crosstalk between adjacent RGC stimulation sites. Photoreceptor integration is achieved using P3HT:PCBM organic photodiodes for cone-associated RGCs and phototransistors for rod-associated RGCs, validated through OghmaNano finite element simulations. The photodiode circuit produces output frequencies of 2.5Hz (dark) to 600Hz (100 W/m2), matching reported RGC response ranges. This architecture eliminates external processing requirements, enabling scalable high-density retinal prostheses design.

High-resolution microscopy techniques are used across research and industry to analyse biological systems, from biomolecules to subcellular organelles, multicellular models and tissues. As multimodal imaging workflows and quantitative analysis of bioimaging data become increasingly widespread, there is a growing need for materials and methods to calibrate imaging systems and evaluate the fidelity of generated image data. Here, we present three-dimensional microscopy phantoms fabricated using two-photon photolithography from transparent resins that exhibit both broadband visible autofluorescence and Raman scattering across the fingerprint and C-H stretching regions. Suitable for analysis using optical profilometry, the phantoms were dimensionally calibrated with SI traceability using a metrological confocal microscope. Immersible in air and common aqueous imaging media, the phantoms are compatible with a wide variety of optical microscopy techniques, including one and two-photon excited fluorescence and coherent Raman scattering microscopy. We employed a forked wedge design to validate image deconvolution results and a stacked lattice phantom to recover image distortion matrices under realistic biological imaging conditions. We demonstrate the impact of correcting chromatic offsets and axial scaling errors for a representative application: analysis of a cell seeded scaffold using confocal laser scanning fluorescence microscopy. These phantoms provide a versatile platform for calibration, quality control and validation of multimodal imaging pipelines and improved quantitative optical microscopy.

Diffusion MRI simulations based on realistic tissue microstructure provide a means to validate biophysical models and optimize acquisition protocols, but their computational cost restricts most studies to domains far smaller than a clinical voxel. The objective of this study was to develop an automated and scalable framework that converts whole-slide histology into diffusion MRI simulations at clinically relevant spatial scales while remaining feasible on standard workstation hardware. We present an end-to-end pipeline integrating two-dimensional whole-slide cell segmentation, mesh generation, and finite element Bloch-Torrey simulation. To enable simulations at large spatial scales without prohibitive memory growth, we introduce a subdomain tiling strategy in which the tissue domain is partitioned into extended subdomains simulated independently under no-flux boundary conditions. Signals are aggregated only from the central regions of each subdomain to minimize boundary artifacts. For an 800 {micro}m x 800 {micro}m histology-based domain, the aggregated signal differed by 0.07% from the corresponding full-domain finite element simulation while reducing wall-clock time from several days to hours and maintaining bounded memory usage independent of global domain size. When applied to a 2016 {micro}m x 2016 {micro}m heterogeneous region approximating the in-plane dimensions of a clinical voxel, the apparent diffusion coefficient obtained from the full domain differed from values computed in smaller dense and sparse subregions, demonstrating the influence of structural heterogeneity at clinically relevant scales on derived diffusion metrics. The proposed framework establishes an automated and memory-stable approach for generating diffusion MRI simulations directly from routine histology.

We find that multi-site temporal control of optogenetic photostimulation in peripheral nerves can enhance firing rates by overcoming the intrinsic limitation of opsin photophysics. The benefits of multi-site optogenetic stimulation were demonstrated with three approaches: (1) in silico modeling, (2) ex vivo in the sciatic nerve, and (3) in vivo in the vagus nerve. An in silico model of multi-site optogenetic stimulation was developed in two Hodgkin and Huxley type neuron models, that supported our hypothesis. The ex vivo sciatic nerve showed an increase in firing frequency that is physiologically relevant for functional control. The technique was then applied in vivo for optogenetic vagus nerve stimulation resulting in significant changes in heart rate compared with standard methods of single-site stimulation. Improving the control of optogenetically induced neural firing will have broad impacts for future developments in optical nerve interfaces and brain-machine interfaces.

Image-based liver Couinaud segmentation is designed to automatically provide the locations of suspicious objects in liver CT/MR images. Once achieved, the physicians will be guided to the target slice and area where the suspicious node is located. However, conventional algorithms trained primarily on healthy liver images often fail to generalize to Hepatocellular Carcinoma (HCC) cases due to pathological structural distortions. In this work, we propose a robust two-stage framework that integrates a 3D Unet with a 3D Anatomical Structure-Guided Graph Convolutional Network (3D GCN). This two-stage strategy effectively isolates the liver volume to eliminate structural noise from neighboring organs, such as the spleen, allowing the framework to focus exclusively on the complex 3D anatomical relationships among the eight segments. To ensure the topological consistency required for global spatial reasoning, we implement a standardized preprocessing pipeline that normalizes liver-only volumes to exactly 50 frames along the z-axis. By combining a lightweight 3D UNet backbone with the 3D GCN for refined boundary reasoning, our model demonstrates superior generalization performance on unseen clinical datasets, achieving a mean Dice score of 0.828 in blind testing. By releasing our code and pretrained weights, we aim to provide the first publicly available deep learning resource for robust Couinaud segmentation.

Altair-dvOPM is an open-access direct-view oblique plane microscope designed for large-field, three-dimensional imaging of cleared and expanded tissue sections. By combining photographic-lens-based detection, externally launched oblique illumination and precision-registered modular baseplates, the system achieves micrometer-scale lateral resolution over a ~5.4 mm field of view without custom objectives or highly specialized alignment procedures. We demonstrate imaging across scales, from subcellular structures in expanded cells to centimeter-scale expanded tissue sections, and provide documentation, CAD files, Zemax models and open-source control software to support replication and extension.

PurposeTo develop a data-efficient deep learning framework for rapid reconstruction of highly accelerated 4D arterial spin labeling (ASL) magnetic resonance angiography (MRA) with robust generalization using extremely limited acquired data, addressing the challenges of prolonged acquisition and reconstruction time. MethodsA simulation-driven, few-shot transfer learning approach was adopted by leveraging publicly available 3D time-of-flight (TOF)-MRA data to generate realistic multi-coil complex-valued pseudo-ASL k-space datasets for large-scale pre-training. A 3D unrolled reconstruction network was trained on this simulated data using a histogram-weighted loss and subsequently extended to 4D using lightweight temporal fusion modules. Fine-tuning was performed using only two experimentally acquired 4D ASL-MRA datasets. The method was evaluated on retrospectively and prospectively undersampled Cartesian 4D ASL-MRA data acquired at 3T and compared with compressed sensing (CS) and locally low-rank (LLR) reconstructions. ResultsThe proposed method achieved superior reconstruction quality compared with CS and LLR, with improved vessel depiction, particularly in distal branches, and enhanced temporal fidelity. Quantitative evaluation demonstrated higher vessel-masked peak signal-to-noise ratio and structural similarity index measure, along with increased error entropy, indicating reduced noise and structured artifacts. The initial pre-trained model already outperformed conventional methods, while additional 4D fine-tuning further improved performance. Robust reconstruction was demonstrated in prospectively undersampled data and multi-slab acquisitions, enabling large-coverage, time-resolved angiography within clinically feasible scan times (4-6 min). ConclusionsSimulation-driven pre-training combined with few-shot fine-tuning enables accurate and rapid reconstruction of highly accelerated 4D ASL-MRA in data-limited settings. The proposed framework provides a practical pathway toward clinically feasible, non-contrast dynamic cerebrovascular imaging.

The retinal nerve fiber layer, composed of axon bundles converging toward the optic nerve, is a key biomarker for diagnosing and monitoring glaucoma and other neurodegenerative diseases. High-resolution en face imaging of individual nerve fiber bundles offers morphological information beyond what conventional optical coherence tomography provides, yet clinical integration remains limited by the lack of automated analysis tools and normative data. Here, we imaged 14 healthy volunteers using time-domain full-field optical coherence tomography and adaptive optics scanning laser ophthalmoscopy, and developed automated pipelines to quantify bundle width, trajectory, tortuosity, and orientation. Bundles were on average 25% wider at shallower retinal depths, width measurements were consistent across imaging modalities, and estimated axon count per bundle decreased significantly with age. Global trajectory analysis revealed systematic deviations of high resolution data from existing mathematical models, particularly in the temporal sector, leading us to propose two refined trajectory models. These normative results provide a foundation for high resolution biomarkers for use in investigations of retinal neurodegeneration.

Tissue phantoms that mimic microvasculature and perfusion are essential for modelling vascular function, guiding interventions, and calibrating imaging systems, which require faithful replication of vascular geometry and flow. Conventional fabrication strategies, including wire-based molding, lithographic micromachining, and additive manufacturing, offer useful capabilities but remain constrained by predefined designs, rectangular channel cross-sections, limited scalability, and high production costs. Reliance on predefined digital vascular models restricts design flexibility and limits the ability to capture the natural variability and complexity of real vascular systems. Here, we present a lithography-free, fractal-generating approach based on a modified Lifted Hele-Shaw Cell (LHSC) technique, in which vascular networks emerge spontaneously via interfacial fluid instabilities. Unlike pre-designed methods, these structures are governed by fluid properties and flow conditions, enabling adaptive, physiologically relevant geometries with smooth Gaussian cross-sections and natural diameter tapering. We demonstrate four phantom designs: a planar vascular tree, an anatomically guided cerebral network, a retinal vascular model, and a conformable curved substrate phantom. Validation using Laser Speckle Contrast Imaging confirms structural fidelity and physiologically relevant flow consistent with Murrays law. This platform uniquely integrates realistic vascular architecture with emergent, fractal driven formation, highlighting its potential as a reproducible and biologically relevant alternative to conventional vascular phantom fabrication. Furthermore, the availability of such realistic in vitro vascular models can reduce reliance on animal experiments and contribute towards more ethical and sustainable preclinical research.